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Abstract 



We suggest a method for integrating sub-families of a family of nonlinear 
SCHRODINGER equations proposed by H.-D. Doebner and G.A. Goldin in 
the 1+1 dimensional case which have exceptional Lie symmetries. Since the 
method of integration involves non-local transformations of dependent and in- 
dependent variables, general solutions obtained include implicitly determined 
functions. By properly specifying one of the arbitrary functions contained in 
these solutions, we obtain broad classes of explicit square integrable solutions. 
The physical significance and some analytical properties of the solutions ob- 
tained are briefly discussed. 



I. INTRODUCTION 

The semi-direct product of the group of diffeomorphisms and the AsELian group of 
smooth functions on M" may be regarded as a generalized symmetry group on M". From the 
representation theory of this group, H.-D. Doebner and G.A. Goldin derived a family 
of nonhnear SCHRODINGER equations on M" [Q-Q], which have been called the Doebner- 
GoLDlN(DG)-equation (For recent progress in the study of these equations see the con- 
tributions in [0.): 

ihdttp = i -—A + v{x) \ ^ + qhDR2[ij]ij + hD'J2cjRAi']i', (i) 

where Rj[ip], j = 1, . . . ,5 are real- valued functionals of the density p := ifjip and the current 
J = Im(V;V?/'), 



V ■ J 

:= A_ 




(2) 



Here the real number D (with the physical dimension of a diffusion constant) labels unitarily 
inequivalent representations of the generalized symmetry group involved in the derivation 
of the nonlinear equations ([l|). It has been interpreted as a quantum number describing 
dissipative quantum systems |^,^ . The real number D' (also with the physical dimension of 
a diffusion constant) describes the magnitude of the real non-linearity and the dimension-less 
constants cj G M are 'model' parameters. 

For the purpose of this paper it is more convenient to use the parameterization that is 
obtained by rewriting equation (|l]) in terms of the real functionals Rj [i/j] only, following the 
notation of 0-0: 

2 5 

F{u, fi) : tdt^p = 1 5^ v^R^^]^ + XI H^A^M^ + Aio^^, yi ^ 0. (3) 
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Particular homogeneous equations of this type have also been considered in the context of 



quantum mechanics by other authors, e.g. |1TT|-[T7 



One of the interesting features of the family of DG-equations (|1]) is its invariance under 
a certain group of transformations [^|| 

N^A^^0) = ',ph^+^+'^^h^-^+^^^ = ItPl e^(7ln|^l+Aarg^)^ 

i.e. if -0 is a solution of F{v,ii), then ■0' = A^(a,7)(V') is a solution of F{u',ii'), where the 
change of parameters under A^(a,7) is 

K = -jl'l + f^U = ^i^l - 7'^2 - -/^l + A/i2 , ^3 = ^> (5) 

f^i = -^/^3 + /i4, /is = ^/^3 - 2/^4 + A/i5, //q = A/iQ. 

Thus, without loss of generality we can restrict our calculations to the particular choice of 
parameters (a particular gauge, see below) 

= -1, i/2 = 0. (6) 

Since the transformations (^) leave the position probability invariant, i.e. p'{x,t) = 
p{x,t), they have been called nonlinear gauge transformations [0,|],|18|. This notion is 
physically motivated by the fact, that in (non-relativistic) quantum mechanics we basi- 
cally measure positions at different times. Furthermore, the transformations have been used 
to construct a consistent notion of observables in a nonlinear quantum theory . 



It turned out that besides such important properties of the DG-equation as homogeneity, 
separability, and Euclidean invariance, which were 'input' by construction, equations (|l]) 
possess a number of other attractive properties. Among them one should emphasize the 
possibility of constructing explicit square integrable solutions, which is important for a 
physical interpretation. In particular, some stationary and non- stationary (GAUSSian and 
traveling wave) solutions have been obtained [^,|6|,p|,pOH2^ . 



The well-known connection between exact solutions of partial differential equations 
(PDEs) and their symmetry properties p5|-p^ as well as the necessity of classifying equa- 



tions (|iD in a unified way, motivated a systematic study of their Lie symmetry in flOl . As a 



result, one has to distinguish nine sub-families (characterized by conditions on the param- 
eters in the chosen gauge) with different maximal Lie symmetry algebras sym'{n). The 
relationship between these sub-families and their symmetries is indicated in Fig.l (using the 
notation of ITI 
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FIGURES 



symo[n) 



/ig— 2/i3 2 

8(^3-1)2 — ~2^1^3' 




Fig. 1: Lie symmetries of the DG-equation. Sub- families are characterized by 
their parameters and arrows indicate the subfamily-structure. The equations 
dealt with in this paper are in bold frames. 



Five of these symmetry algebras are finite dimensional, symo(n) - sym^^in). Among 
them we find the direct sum of the (centrally extended) Schrodinger algebra and the 
real numbers (due to real homogeneity of the equations). These equations thus fit into the 
classes of Schrodinger invariant nonlinear evolution equations determined in [p8|-|31|. 

The four remaining symmetry algebras are infinite dimensional. sym\{n) and sym\{n) 
contain in addition to the elements of symi{n) infinite dimensional algebras h and c, that 
depend on a pair of (real) solutions of a linear forward and backward heat equations and 
a (complex) solution of a linear Schrodinger equation, respectively. Actually these sym- 
metries correspond to linearizations of these sub-families, the first to a pair of forward and 
backward heat equations, the latter to a Schrodinger equation On the contrary, 

the symmetry algebras sym'^in) and symQ{n) contain an infinite dimensional algebra a that 
depends only on one real- valued function. As a consequence, there is no local transformation 
(i.e. a transformation that does not involve integrals or derivatives of the dependent vari- 
ables) linearizing the corresponding DG-equations. Nevertheless, these equations as well 
as the one admitting the finite dimensional symmetry algebra sym^{n) are shown in the 
present paper to be integrable by a non-local transformation of dependent and independent 
variables in the case of one spatial variable (n = 1). Thus, all DG-equations with excep- 
tional symmetries (bottom row of Fig. 1) sym\, sym\, sym^, sym2, syrriQ are integrable, i.e. 
they can be reduced to an equation which is either linear or integrable by quadratures. 
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The principal object of study in the present paper are DG-equation in 1+1 dimensions 
with parameters 



Z/i = -1, U2 = 0, /ii = /i2 = ^4 = ^5 = 0, (7) 
Vl = -1, Z/2 = 0, /i2 = /is = 0, /i3 = 2, /i4 = -/ii 7^ 0, (8) 

i.e. the following coupled two-dimensional PDEs: 



=lif^i- ^)Im ^ + 2 (^Im^^ ' - /.Jm (^^^ ' + i^oVix) i ^. (10) 



One of these DG-equations is contained in the so-called Ehrenfest sub-family ful- 
filling the second Ehrenfest relation: equation (P) with /i3 = 1, i.e. the DG-equation with 
maximal Lie symmetry sym^{n). This equation is furthermore the only Schrodinger and 
therefore Galilei invariant equation among (^) and ([10|) . Nevertheless, in the free case 
[V = 0) all of these DG-equations admit traveling (solitary) wave solutions with arbitrary 
shape 0. These solutions are rediscovered as a particular case of the general solutions in 
this paper. 

Using a polar decomposition 

'ilj{x,t) = exp {r{x,t) + is{x,t)) (11) 

we rewrite the above equations in the following way: 

^^in + s.. + 2r.s = 0, 

^^(n + s^. + 2r^s = 0, 
[ st + /iis^x + 2s^ = -jJoV. 

The paper is organized as follows. In the section |I| we integrate the free equations, i.e. 
equations (0), 0T3| ) with a vanishing potential (V = 0). In order to integrate Fi we have to 
distinguish between the cases /xs 7^ 1 and /X3 = 1, the latter corresponding to the subfamily 
with the larger Lie symmetry algebra sym^{n) D sym2{n). 

Section |T| contains some remarks on the integration of the equations with potential and 
two particular examples where the integration is carried out. The methods of integration of 
Fi and F2 in sections y and |T1| yield their general solutions containing implicitly determined 
function. Consequently, these solutions are, generally speaking, implicit. Therefore, in 
section |IV| we give some explicit solutions for the free equation as well as for linear and 
quadratic potentials by specifying one of the arbitrary functions of the general solutions 
obtained in the preceding sections. 
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II. INTEGRATION OF FREE DG EQUATIONS 



Putting in (|T3p V = we obtain the following PDEs: 



1"t ~l~ ^xx ~l~ "^^x^x 
St + ^3^ 

~l~ ^XX ~l~ '^^X'^X 

St + fJ-lSxx + '2sl 



0, 
0, 

0, 
0. 



(14) 
(15) 



Henceforth we suppose that in (^5|) /ii 7^ 0, since otherwise system (|1^) is a particular case 
of (|14D with /is = 2. 



A. Integration of the family Fi 



First, we turn to the integration of the system of nonlinear PDEs ([T^. As this system 
admits only a finite-dimensional Lie symmetry group, there is no local transformation which 
linearizes it. So the only possibility to transform the system in question into an integrable 
form is to utilize a non-local transformation of dependent and independent variables. The 
choice of a desired change of variables is implied by the form of the second equation of the 
system (p!4D; it is nothing but the one-dimensional HAMILTON- J ACOBI equation, which is 
known to be linearizable by the following contact transformation (called in the literature 
the EULER- Ampere transformation): 

Zo = t, Zi = Sx, U = XSx- S, = -St, Uz^ = X. (16) 

Let us recall that a transformation (x, t, s(x, t)) ^ {zq, z\, ^(-ZO) -^i)), 

fcXs^^ ^) "5) Sx-i St), 

u = g[x,t,s,Sx,St), (17) 

h(^(^X, t, S, Sx, "St), 

where a = 0, 1, is called contact, if it preserves the first-order tangency condition 

ds — Sxdx — Stdt = =^ du — Uz^dzQ — Uz^dzi = 0. 

The above condition ensures that the functions u^^ , Uz^ determined by the last two formulae 
from (|l^) are really derivatives of a function u determined by the third formula with respect 
to Zq and Zi determined by the first and second formula, correspondingly. 

It is readily seen that formulae (|16D determine a contact transformation preserving the 
tangency condition. But before applying the Euler- Ampere transformation to the system 
under study we must ensure its invertibility, as we may loose some solutions otherwise. It is 
known that transformation (plGf ) is invertible in a domain where Sxx 7^ 0. Consequently, we 
have to consider the cases Sxx ^ and s^x = separately. 
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Case 1. Sxx ^ 



Let us apply the transformation ( ]TE| ) to the system (|Tip having prolonged it to the 
second derivatives 



Six ^ttSxx ^tx 1 



' ^ ■ (18) 

^ 3)X ^ XX ^ XX 

As a result, we get 

^i^o = ^^3%, ^^20 H ^'^i = • (19) 



Now, the first equation becomes linear and is easily integrated to give the following 
expression for ^(^0,-21): 

u = fX3Zozl + f{zi), fEC\R,R). (20) 

Inserting the result into the second equation of the system (^) we get a first-order 
linear PDE with non-constant coefficients, 

r I 2(l-/i3)zi ^ 1 

When integrating the above equation we have to distinguish two sub-cases = I and 
/i3 7^ 1. Let us recall that the DG-equation with parameters under /ia = 1 satisfies 
the Ehrenfest relation and, what is more, admits an additional symmetry operator 
(see Fig.l). 

Sub-case 1.1. /ia = 1 



In this case equation (21) takes the form 

1 



2^0 + f'izi) 
and its general solution reads 

r{zo, z,) = -\ ln(r (zi) + 2z,) + g{z^), g e C\R, R). (22) 

To rewrite the result obtained in the initial variables {x,t,r{x,t), s{x,t)) we have to 
invert the transformation (|16|). From the second and third relations it follows that 

s = XZi — u = XZi — tz\ — f\zi). 



To determine the function zi = zi{x,t) we make use of the last relation from (|T6|). 
Substituting into it the formula ([20|) with /^a = 1 yields 



X = u^, = 2zoZi + f'{zi), 



hence 

2tzi- x + f'{zi) = 0. (23) 

The above relation determines the function zi{x,t) in an imphcit way. Since s^x = 
2t + f"{zi) ^ 0, we can always solve (p3|) (at least locally) with respect to Zi thus 
getting an explicit form of the function zi. 

Summarizing the results we conclude that the general solution of the system of PDEs 
(|14D with ^0, /i3 = 1 is of the form 



r{x,t) = -iln(/"(^i) + 2t) + (7(zi) 
s(x,t) = —tz\ + xzi — f{zi)^ 



(24) 



where f,g & C^(M,M) are arbitrary functions and zi = zi{x,t) is determined by the 
relation (ESl). 



Sub-case 1.2. /xs 7^ 1 

Using in (21) the transformation 



r{x,t) = r{x,t) - ^ Inzi 

we get 

. 2(l-/i3)zi . _ 

'"20 o \ till 

2/13^0 + / (^l) 

From the theory of the first-order PDEs it is well-known (see, e.g. [^) that a general 
solution of the above equation has the form f = giyj^zi^^zx)), where u){zi^,z\) is an 
integral of the Euler-Lagrange system 

dzQ dz\ 



I 2(1-/13)2:1 ■ 
2^l3zo+f"izl) 

The above system is rewritten as a linear first-order ordinary differential equation for 
a function Zq = Zq{zi), 

"^'^ ^ '2^3^o + r(^i) 



dzi 2^1(1-/23) 
the general solution of which can be represented in the form 
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c = 2(^3 - i)zozr" + / c^r(c)< 



with an arbitrary constant C. Hence we conclude that the general solution of equation 
( PID is given by the following formula: 

r(x,t) = ^^^^^— ^ln^i + (7f2(/i3-l)^o^r"+ j (^/"(O^^Cj , (25) 

where f,g& C^(M,M) are arbitrary functions. 

Returning to the initial variables x, t, r, s we obtain the general solution of DG- 
equation (U) for the case Sxx ^0, 7^ 1 



s(x,t) = -fistzf + XZi - f{zi), 



(26) 



where f,g E C^(M, M) are arbitrary functions and the zi = zi{x,t) is determined 
implicitly 

2fi3tz, -x + f'{zi) = 0. (27) 

Case 2. s^x = 

In this case s{x,t) = a{t)x + f3{t) with arbitrary smooth functions a{t), j3{t). Substi- 
tuting this expression into the second equation from ([I^) we arrive at the relations: 

a'{t)x + (3'{t) + /i3a^(t) = 0, n + 2a(t)r^. = 0. 

An integration of these equations gives rise to the following expressions for r and s: 

r{x,t) = f{x-2Cit), 

(28) 

s{x,t) = Cix - ^isCft + C2, 

where Ci,C2 are arbitrary real constants. Thus we have rediscovered the traveling 
(solitary) wave solutions with arbitrary shape as a particular case of the general 
solution. 

We have established that any smooth solution is contained (at least locally) in one of the 
classes given by equations (j2^), (0), and (pS]). Summarizing we arrive at the conclusion 
that the general solution of the free DG-equation (P) splits into two inequivalent classes: 
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1. /X3 = 1 

'ilj{x,t) = f{x - 2Cit) expf^i{Cix - C^t + C2)Y (29) 
tP{x,t) = (/"(^i) + 2ty'giz^) exp{-t[tzl - xz, + /(^i)) }, (30) 

where /, g are arbitrary sufficiently smooth functions, Ci, C2 are arbitrary real param- 
eters, and Zi = Zi{x,t) is determined implicitly by formula (p^); 

2. /Z3^1 

^(x,t) = f{x-2Cit)expJ^i{Cix- i23Cft + C2)y (31) 



ij{x,t) = z^^'-^-'^g(2{f,,-i)tzr' + j (—nodc) 

X exp|-z(^/i3t2;2 - xzi + /{zi)^, 



(32) 



where f,g are arbitrary sufficiently smooth functions, Ci, C2 are arbitrary real pa- 
rameters, and zi = zi{x,t) is determined implicitly by formula (p^). 

Although formulae (|3l|) and (|3^ ) give the general solution of the corresponding DG- 
equation for all /is 7^ 1, the case /is = deserves a special consideration, as the system of 
PDEs ([1^ with /xs = is easily integrated without applying the contact transformation 
(p!6|), (0). The second equation of (|T^ yields that s does not depend on time, so 

s{x,t) = f{x), feC\R,R). 

Inserting this into the ffist equation we get a ffist order PDE for r, 

n + 2f'{x)r, + fix) = 0. (33) 

The case s{x,t) = const leads to a traveling wave solution (^) with /ts = 0, Ci = 0; if 
f'{x) ^ 0, then the general solution reads 

r{x,t)=g(^2t- f J^yj-l^^f'^^)^ (34) 

where g G C^(M, M) is again an arbitrary function. 

Consequently, the general solution of the DG-equation (P) with /is = is either given 
by a traveling wave solution ( |3TD with /is = 0, Ci = 0, or by 

^(x, t) = [fix)) '^'g {21 - exp{z/(a:)}. (35) 
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B. Integration of the family F2 



Let us turn to the integration of the DG-equation (|15|) . First, we note that the second 
equation is the potential Burgers equation, which is is hnearized by the logarithmic substi- 
tution. Furthermore, we reduce the order of spatial derivatives in the first equation by a 
linear transformation of the dependent variables. Thus, the transformation 

t) = — ;Uir(x, t) + t), 'u(x,t) = exp( — (36) 

reduces system ([TsD to the form 

UVt + UlUxV^ = 0, Ut+ /ilMxx = 0. (37) 

The second equation may be taken as the integrability condition of the vector-field {u, —fiiUx) 
on space-time, dtu = dx{—f^iUx), so that it is the gradient of a smooth function (f{x,t), 

Lpt = -fiiu^, (fx = u. (38) 
With this remark the first equation of the system ( |57D is rewritten to be 

ifxVt - (ptVx = 

and is easily integrated v{x,t) = f{{p{x,t)), where / is arbitrary, sufficiently smooth function. 
Solving ( PSP with respect to ip we get 

ip{x,t)= ! u{^,t)d^ - Hi [ Ux{0,T)dT + C, (39) 
Jo Jo 

where C is an arbitrary constant. 

Returning to the initial variables {r{x,t), s{x,t)) we get the general solution of the DG- 
equation (p!0| ) 

^ix,t) = {u{x,t))-^f(^J^ u{^,t)d^~^iiJ^Ux{0,T)dT^exp\^'-^\nu{x,t)y (40) 

where u{x,t) is an arbitrary solution of the heat equation 

ut + fJ'iUxx = 0, (41) 

and / is an arbitrary smooth function. 

Finally, we note that the traveling wave solutions of are reobtained using the particular 
solution 



u(x, t) = exp <( {x — vt) 

Hi 



of the heat equation (^Tj). 
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III. DG EQUATION WITH NON- VANISHING POTENTIAL 

Surprisingly enough, DG-equations (|I^), (|T^) are integrated in quadratures even in the 
case when V{x) ^ (i.e. in the presence of a non- vanishing potential). Unfortunately, for 
the family Fi the corresponding formulae are implicit and cumbersome. That is why we 
restrict ourselves to considering in detail system (O) with an additional constraint /xs = 1 



i.e. the Ehrenfest subfamily of (12) is studied); this system was also considered in |23 



A. Integration of the family Fi 

Choosing in (|12|) /ia = 1, we obtain the following system of PDEs: 

n + 2r^s^ + s^^ = 0, st + si + ^ioV{x) = 0. (42) 



To linearize this system we make use of the following trick: instead of system of PDEs (^21) 
one of its differential consequences is considered 

rt + 2r^s^ + s^^ = 0, Stx + 2s^s^^ + fioV'{x) = 0. (43) 

Substituting in (|43|) 

R{x,t) = r{x,t), S{x,t) = s^{x,t) (44) 
we arrive at the system of first-order PDEs 

Rt + 2SR^ + S^ = 0, St + 2SS^, + fioV'{x) = 0. (45) 



Thus, using the substitution (^J) enables us to reduce the order of system of PDEs under 
study. Next, we apply to this system the hodograph transformation (x, t, R{x, t), S{x, t)) i-^ 
{zo,zi,v{zo,zi),u{zo,zi)) 



t = U, X = Zi, S = Zq, R = V, 



R. = ^ Jl =y _ Si = — S = ^^^"^ 



This hodograph transformation is defined in an arbitrary domain where St ^ 0. So again 
we have to distinguish two cases, St ^ and St = 0. 

Case 1. St^O 

Performing in (EH) the change of variables (EB) we get 



u 

fioV'{zi)u^,^ - 2zou^, + 1 = 0, ^0 ^'(2^1)^20 - '^^oVzi + — = 0. 

A further change of variables 

u{zo,Zi) = u{zo,Zi), v{zo,Zi) = v{zo,Zi) + llnu^^{zo,Zi) (47) 
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transforms the system to 

/ioV" - 2zqU^^ + 1 = 0, ^^oV'{zl)v^Q - 2zoV^^ = 0. (48) 

Thus, combining local and non-local transformations of the dependent and independent 
variables we linearized and decoupled the differential consequence of (|42D . Integrating 
these equations yields 



uizo,z,) = i P (^zl + fio(viz^)-V{0')y' dC + f(^z^,+fioViz,)y, 
v{zo,zi) = gl^z'^ + fioVizi)^ 



(49) 



where f,g& C^(R,M) are arbitrary functions. 

Returning to the variables (^x,t, R{x,t), S{x,t)j we obtain the general solution of 
system (^) in an implicit form 

^ = i (^'(^' + ^0 {Vix) - ViO)) '"di + f {S\x, t) + ^oV{x)) , (50) 
R{x, t) = \\n St{x, t)+g (^S\x, t) + fioV{x)^ . (51) 

To rewrite these equations in the initial dependent variables {r{x,t), s{x,t)) we have 
to invert the transformation (0). As a result, we get 



X 



r{x,t) = R{x,t), s{x,t) = S{T,t)dT + Lp{t), (52) 





where Lp G C^(R,M) is an arbitrary function. The ambiguity arising is connected to 
the fact that we are not solving the initial system (^2]) but its differential consequence 
(^3]). The 'extra ' function (p{t) is used to choose from the set of solutions of system 
( ^51) those ones which satisfy (^). Indeed, substituting formulae ( ^2]) into (^) and 
taking into account that the functions R{x,t), S{x,t) satisfy system of PDEs (^) we 
arrive at the following ordinary differential equation for a function f{t): 

so 

^{t) = - [ S\0, T)dT - fioV{0)t + C, 

Jo 

where C is an arbitrary real constant. 

Summing up, we conclude that the general solution of the initial DG-equation reads 
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r{x,t) = -^\nSt{x,t) + g(^S'^{x,t) + fioV{x)^, 
s{x,t) = / S{T,t)dr- / S'^{0,T)dT- fioViO)t + C, 



where S{x, t) is a smooth function determined imphcitly by (^) and /, g are arbitrary 
sufficiently smooth functions. 

Case 2. St = 

With this condition the system of PDEs (H3p is easily integrated to yield 



/3 

S{x,t) = a/Ci - fioV{x), 



where g is an arbitrary sufficiently smooth function and Ci is an arbitrary real constant. 
Rewriting the above expressions in the initial variables {r{x,t), s{x,t)) we get 

d^ 



r{x,t) = -\\n{C,-fioV{x))+g(2t- j -j= 



s{x, t) 



r VCi - MOd^ - Cit + C2 

Jo 



(54) 



where C2 is an arbitrary constant. 

Thus, we have established that the general solution of DG-equation (P) with /is = 1 
splits into the following two classes: 



1 

i{j{x,t) = (^St{x,t)y g(^S\x,t) + fioVix)') 

X exp I* (^J^ S{C, t)dC - S^O, T)dT - fXoV{0)t + C 



(55) 



where g G C^(]R, M), C G M, and S{x,t) is determined implicitly by (^Ol). 
2. St = Q 

di 



xexp|^(^^ ^Ci- ^ioV{i)d^-Cit + C^'^ 



(56) 



where g G C^{^, M) and Ci, C2 G M. 
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B. Integration of the family F2 



Integrating system (|T3D with potentials is similar to integrating the free system {V{x) = 
0) in section |II B| . Using the change of variables (|36|) for (|13D we arrive at the following 
system of PDEs for new functions u{x,t) and v{x,t): 

ut + ^^lU^x + — = 0, uvt + ^llU^Vx + ^qV{x)u = 0. (57) 
/ii 

Now, given a potential V{x) and an arbitrary solution u{x, t) of the first equation of this 
system, one can construct a general solution v{x^ t) of the second equation which leads to a 
general solution of the initial system: 

1 

tp{x, t) = (^u{x, t) j ^ exp (^—fi^^v{x, t) + In u{x, t)j . (58) 



IV. EXPLICIT SOLUTIONS 

As mentioned before some of the solutions of DG-equation obtained are local in a sense 
that they are not determined on the whole plane M^. But for physical applications one needs 
global solutions, and what is more, they should be square integrable, i.e. the integral 

/oo 
ilj{x,t)ip{x,t)dx (59) 
-00 

is to be finite. If it is, the quantity p{x,t) = ^ip{x,t)'i{j{x,t) > is treated as a probability 
density of a distribution of the wave function i(j in space at a given time. 



A. Explicit solutions of the family Fi 

Evidently, the traveling wave solutions of DG-equation (|^) given by (^9D, (|31|) are defined 
on the whole plane and, consequently, are global. To ensure square integrability of these 
solutions one has to restrict the choice of the arbitrary function / to square integrable ones, 

/oo 
/^(r)(ir < 00. 
-00 

Thus, the traveling wave solutions are square integrable provided / is. 

Solutions (|30|), ( ^2]) are, generally speaking, local, since the function zi{x,t) contained in 
these solutions is determined implicitly by formulae (^31) and (^7\), correspondingly, and the 
existence of solution is only guaranteed locally by the implicit function theorem. In order 
to obtain explicit expressions for global and strictly local solutions we consider solutions of 
(pop with quadratic and cubic functions /, respectively. 

For quadratic functions /, 

f{zi) = fisazl , 
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the implicit equation ( pT)) (resp. (0)) for zi can be solved globally and we get zi{x,t) ~ 
2^3 (t-a) • Thus, we arrive at the following class of explicit solutions of the DG-equation @ 
containing an arbitrary smooth function g: 

1. /i3 = 1 

\ f ■ 2 

tXJ \ I LtXj 



2. /13 7^ 1 



^ix,0 = (2r-9(^(^jexp|j^^ (60) 



• 2 

IX 



ijj{x,t) = {t- Q;)2(i-M3)a:2(A'3-i)^ {t - "3-1 X '^3"! exp <^ \ . (61) 

These solutions are square integrable, provided g is, and are well defined on the whole plane 
with a possible exception of the line t = a, where they converge to a 5-function. In 
particular for g{z) = exp(— 2^) solutions (|60D coincide with the GAUSSian wave solutions of 

El- 

Cubic functions /, 

give rise to strictly local solutions of the DG-equation. Indeed, inserting them into ( p3D we 
obtain a quadratic equation with respect to zi 

zl + 2/i3tzi - X = 0. (62) 

It has real solutions in the case x + > only, i.e. z\{x,i) is not defined inside the 
parabola x + = 0. Solving (|6^) yields z\{x,i) = — t ± y^x+Tip^; according to the 
general solutions (^) and ( P^ we have to choose the positive sign since in (^) It + j"{z\) 
and in ( p2D z^ have to be positive. Hence, we arrive at the following class of strictly local 
explicit solutions of the family Fx (|^) containing an arbitrary smooth function g: 

1. /i3 = 1 



3 



?/'(x,t) = (x + t^Y'^g [-t + + exp |-*(^^^ + - ^(x + t^)5 ) J> , (63) 



2. /i3 7^1 



^(x,t) = (y^T#-/i3)^'''^-%(^(v/^T#-/i3^ 

X (^(Ai3 - l)t + v/^T#) ) (64) 
X exp 1^ |^2(x + ii\t^){-ii^ + sJx + Yif^) + /i3tx^ I . 
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The domain of definition of these solutions is the set {{x, t) : x + fi^t"^ > 0}. Furthermore, as 
the function a/x + js^t"^ is not defined for x at — oo at any given time t, \t\ < oo, the solution 
(|63|) is not square integrable. 

In this context let us remark that in general solutions (^) are square integrable at a 
given time t provided 

• the (possibly infinite) limits a± = lim zi{x,t) exist and 

z— >dzoo 

• g is square integrable on the interval [a_, a+]. 

This statement follows from a change of the integration variable x Zi{x,t). 

Before turning to DG-equations with non-vanishing potentials we examine solution (|35D 
of the particular case /is = of the family Fi. If the first derivative of / has no zeros, then 
the solution given by (^5]) is certainly global. Again, a change of the integration variable 
shows that the solution is square integrable, provided that 

r dr 

• the (possibly infinite) limits a± = lim / exist and 

• g is square integrable on the interval [a_,a+]. 

In case of non-vanishing potentials we concentrate on the following specific potentials: 

1. the linear potential 

V{x) = —x, aeR] (65) 

2. the harmonic oscillator potential 

l/(x) = — x^ aeM; (66) 

3. the anti-harmonic oscillator potential 

V{x) = x\ aeR; (67) 

First we consider Ehrenfest case /is = 1, the integration of which has been studied in 
detail in section |IIl7\| . 

1. For linear potentials (|65|) the implicit equation (|50D reads 

t = --S{x,t) + f{S^{x,t) +ax). (68) 
a 

If we choose / = 0, then S{x,t) = —at. Thus, we get a class of explicit solutions from 
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These solutions are defined on the whole plane and square integrable, provided g 
is. 



Another class of explicit solutions is obtained directly by means of formula (|56|): 

t) = (Ci — ax)^ig (^at + v^Ci — ax^ exp ~ ^^)^ ~ ~^ , 

(70) 

where g is an arbitrary twice continuously differentiable function, Ci, C2 are arbitrary 
parameters. These solutions are defined on the half-plane {{x,t) : x < ^}. 

Analogously we construct explicit solutions for the ( ant i-) harmonic oscillator potentials. We 
give these without derivation. (/ is an arbitrary sufficiently smooth function, Ci, C2 are 
arbitrary parameters.) 



2. harmonic oscillator potential ( |66D 

ip{x, t) = x~2 f (x~^ sin 2at) exp |^(^^^^ 2at + Ci^ | , (71) 

ipix, t) = {Cf - a^x^)--^f (^^Jcl-a^x"^ sin 2at - ax cos 2at^ (72) 

X exp V ( o V '-^i ~ '^^^^ + ^ arcsin — - Cit + C2 
^^ \2 ^ 2(1 



3. anti-harmonic oscillator potential (|67|): 

?/^(x, t) = x'^f [x^^ sinh 2at) exp |^ (^^^ coth 2at + Ci j | , (73) 

tp{x,t)= (Ci + a^xYKf [[ax + y^d + a^x^^ e""'^*) (74) 

X exp a/Ci + a^x^ -|- — In ax + a/Ci -|- a^x^ 
[ V2 2a 



Cit + 



As mentioned in section |111 A| , general solutions of the family Fi (^) with fi^ ^ 1 are 



given by cumbersome implicit formulae. But with the particular choice of the potentials 
above it has explicit solutions containing one arbitrary function: 

I. fi3 = 

1. linear potential (|H3|): 

ip{x, t) = g{x + at^) exp{— zatx} ; (75) 
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2. harmonic oscillator potential ([66|): 

3. anti-harmonic oscillator potential (|67|) : 

ip{x,t) = e""'*'^ (^xe"^"'*') exp{m^te^}; 

11. = > 

1. linear potential 



tp{x, t) = g(x + at^) exp < —i ^atx H g"^^ ~ 



2. harmonic oscillator potential (|66|): 

ip{x,t) = x~^g ^x"^^ sin 2aAt j exp j^^^a;^ cot 2aAt + | ; 

3. anti-harmonic oscillator potential (|67|) : 

'(/'(x, t) = x~2g ^x~^^ sinh 2aAtj exp |^(^^'''^ "^^^^ 2aAt + ^) | 

III. /is = -A^ < 

1. linear potential (p5|): 



t) = (^(x + at^) exp I —2 (^atx '^^^ ~ ^) j" ' 

2. harmonic oscillator potential (pB|): 

t) = x~2 ^x'^^ cosh2aAtj exp |^(^^^^ tanh2aAt + | 

3. anti-harmonic oscillator potential (|67|) : 

t) = x~ 2 ^x'^^ cos 2aXt^ exp |^ (^^^ 2aAt + C 

In all these cases / G ^^(M, M) and C g M. 
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B. Explicit solutions of the family F2 



Clearly, if the function u{x,t) is a global solution of the heat equation (^T]), then the 
formula (|40D gives a global solution of DG-equation (|10]). And what is more, it is square 
integrable provided 

• the (possibly infinite) limits a± = lim I / u{T,t)dT — /ii / Ux{0,T)dT ) exist and 

±~ \Jo Jo ) 

• f is square integrable on the interval [a_ , a+] . 

In order to construct explicit solutions of the family F2 (|13D for linear and quadratic 
potentials (|53|)-([H7|) we have to solve equations (^T]). After some tedious calculations we 
obtain the following solutions: 

1. linear potential (|65D: 

C 2ia 1 

tpi^x, t) = f\x + at^) exp < -iatx —t^ + iC \ ; (84) 

2. harmonic oscillator potential (|66D: 

ipixjt) = (cos4at)^2/ (x^^ cos4af) 

X exp I —i (^^3;^ tan Aat + ^ In cos Aat — C j | ; ^^^^ 



3. anti-harmonic oscillator potential (|67|): 

ip{x,t) = (cosh4at)~2j (^a;~^ cosh4at) 



X exp |^(^^^^ tan4at — ^ In cosh4at + icj | . 



(86) 



Here / is an arbitrary twice continuously differentiable function and C is an arbitrary con- 
stant. 



V. CONCLUSION 

As mentioned above Lie symmetries of DG-equations considered are not extensive 
enough to provide their linearizability by means of local transformations. PDEs (P), (p!oD 



prove to be integrable because of infinite non-local symmetries admitted. Take, as an exam- 
ple, system (^21). It has been decoupled into a system of two linear first-order PDEs (|48D by 
means of non-local transformations of dependent and independent variables (0), (^Bj), {^7\). 
It is well-known (see, e.g. |^^) that any linear first-order PDE admits an infinite parameter 
Lie transformation group. Consequently, system ( ^H] ) possesses an infinite local symmetry. 
But after being rewritten in the initial variables {x,t,r{x,t), s{x,t)) it becomes non-local 
and can not be found by using the infinitesimal Lie algorithm. 
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In the case involved an existence of non-local symmetry was indicated by a change of 
local symmetry of the DG-equation when the parameters were specified to be (0), (|]). These 
additional local symmetries form the top of the 'iceberg', the main part of which consists of 
non-local symmetries enabling us to integrate the corresponding DG-equations. 



Since we have the formulae for general solutions of systems of PDEs (^), (|iy), it is not 
but natural to apply these to analyze the initial value problem for these systems, which is 
important for a physical interpretation of the equations. For example, using formula ( |1D| ) it 
is not difficult to prove that the initial value problem 

iift = Ufil- i)lm 




i/j^x, 0) = ro(x) exp{zso(x)}, 

where ro,so G C°°(]R, M) are arbitrary functions such that /iiSo(x) 7^ —00, has a unique 
solution given by the formula (^OD, where u{x,t) is a solution of the initial value problem 
for the heat equation 

Ut + fJ'iUxx = 0, u{x, 0) = exp I —So{x] 
and the function f{y) reads 

/(y) = /?(%)) exp 

Here h{y) is determined implicitly by the relation 

rKy) r 2 "1 
y = exp <^ —So{r) > dr. 



L A'l 

But for the DG-equation (P) an analysis of the initial value problem is complicated due to 
the complex structure of its general solution. 

The method of integration of DG-equations developed in the present paper for the case 
of one space variable can be extended to a physically more interesting case of three spatial 
dimensions. A principal idea of such an extension is a utilization of generalized Euler- 
Ampere transformations of the space (t,x,u(t,x),Ut(t,x),gradu{t,x)) suggested in P2| . 
The above transformations were used to study compatibility and to construct a general 
solution of the four- dimensional nonlinear D'ALEMBERT-eikonal system. This problem is 
under investigation now and will be a topic of our future publications. 
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